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Preface 


Although  the  phenomenon  of  laser  induced  gas  break- 
down has  been  extensively  studied  for  several  years,  lit- 
tle progress  has  been  made  in  developing  a simple  model 
which  can  quickly  and  accurately  predict  breakdown  at  low 
values  of  incident  flux.  The  investigations  described  in 
this  thesis  represent  an  attempt  to  develop  such  a model. 
While  success  was  limited,  it  is  hoped  that  the  informa- 
tion presented  here  might  prove  to  be  a valuable  starting 
point  for  anyone  desiring  to  continue  investigation  of 
this  problem. 

I would  like  to  extend  my  sincere  appreciation  to  my 
advisor,  Capt.  P.  E.  Nielsen,  who  suggested  this  problem 
and  provided  the  guidance  needed  to  successfully  complete 
this  project.  I would  also  like  to  thank  Mrs.  Mildred  C. 
Kelley  for  her  aid  in  the  proofreading  and  preparation  of 
the  final  report. 

Finally,  my  deepest  appreciation  is  extended  to  my 
wife,  Mary,  who  was  a source  of  encouragement  and  inspi- 
ration not  only  during  this  project,  but  throughout  my 
eighteen  months  of  study  at  AFIT. 
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Abstract 


An  attempt  is  made  to  develop  a simple  yet  accu- 
rate gas  breakdown  model  which  can  be  easily  coupled  to 
the  hydrodynamic  equations  governing  fluid  flow  in  laser- 
target  interactions.  The  accuracy  of  three  relatively 
simple  models  is  investigated.  Each  is  compared  with  the 
more  accurate  and  complex  quantum  kinetic  model,  to  deter- 
mine the  conditions  under  which  it  maintains  reasonable 
accuracy.  A gas  consisting  of  a single  monatomic  species 
is  assumed  and  attention  is  restricted  to  the  early  por- 
tion of  the  electron  cascade.  A temperature  model  is 
found  to  agree  reasonably  well  with  the  quantum  kinetic 
model  at  values  of  incident  laser  flux  greater  than  5 x 
107  W/cm  . A diffusion  model  is  found  to  yield  similar 
results.  A two-temperature  model,  which  is  derived  in 
an  attempt  to  extend  the  range  of  the  temperature  model 
to  lower  values  of ' incident  flux,  is  found  to  be  invalid. 


vii 


I.  Introduction 


The  breakdown  threshold  of  the  atmosphere  imposes  an 
upper  limit  on  the  amount  of  energy  which  can  be  propagated 
through  the  atmosphere  by  a laser.  Below  this  rather  well 
pronounced  threshold  intensity,  the  atmosphere  is  relative- 
ly transparent  to  the  laser  beam.  However,  once  the  thres- 
hold is  exceeded  cascade  ionization  ensues.  Naturally 
occurring  electrons  are  rapidly  heated  by  the  laser  radi- 
ation. Those  attaining  sufficient  energy  ionize  neutral 
species  and  create  more  electrons.  The  net  result  is  the 
formation  of  a highly  ionized  plasma  which  severely  atten- 
uates the  beam.  Similarly,  vapor  emerging  from  a target 
irradiated  by  a laser  can  break  down  and  shield  the  tar- 
get from  further  deposition  of  laser  radiation. 

Under  the  sustained  action  of  a laser,  several  inte- 
resting post-breakdown  effects  will  emerge.  The  newly 
formed  plasma  continues  to  absorb  almost  all  the  energy 
incident  on  it,  heating  the  gas  in  the  absorption  region 
to  very  high  temperatures.  As  this  relatively  small  vol- 
ume of  gas  begins  to  expand,  the  surrounding  air  is  heated 
and  becomes  ionized.  This  causes  a new  layer  of  gas  capa- 
ble of  absorbing  energy  to  be  formed  within  the  beam's 
column,  and  a "laser  absorption  wave"  will  propagate  to- 
wards the  laser, decoupling  it  from  the  target.  At  very 
high  laser  intensities,  the  hot  gases  will  expand  so  ra- 
pidly that  a shock  wave,  capable  of  interacting  with  the 
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target, will  be  formed. 

If  these  hydrodynamic  phenomena  associated  with  laser- 
target  interaction  are  to  be  properly  analyzed,  an  accurate 
prediction  of  the  breakdown  threshold  is  required.  Present 
calculations  of  breakdown  phenomena(the  quantum  kinetic 
model)  use  lengthy  computer  programs  which,  even  when  sim- 
plified, require  an  average  of  30  seconds  for  each  run. 
These  calculations  must  be  repeated  at  least  300  times  for 
each  centimeter  the  shock  wave  travels.  Consequently, 
there  exists  a need  for  a simple  yet  accurate  breakdown 
model  which  could  be  easily  coupled  to  the  equations  gov- 
erning fluid  flow  in  a laser-target  interaction. 

In  an  effort  to  find  such  a model,  three  relatively 
simple  breakdown  models  were  investigated.  The  accuracy 
of  each  was  determined  through  a comparison  with  the  more 
accurate  quantum  kinetic  model.  Primary  emphasis  was 
placed  on  comparisons  involving  a temperature  model.  This 
model  assumes  a maxwellian  distribution  of  free  electrons. 
The  second  model  compared,  a diffusion  model,  assumes  that 
electron  heating  can  be  represented  as  a diffusion  of  free 
electrons  along  the  energy  axis.  Finally,  a two-temper- 
ature  model  which  assumes  a segmented  electron  distribu- 
tion function  was  developed  and  its  results  compared  to 
those  of  the  quantum  kinetic  model. 

A brief  review  of  the  theory  involved  in  laser  induced 
gas  breakdown  is  presented  in  Chapter  II.  Chapter  III 
provides  a discription  of  both  the  comparisons  made  and 
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the  models  compared.  Results  are  presented  and  discussed 
in  Chapter  IV  followed  by  conclusions  and  recommendations 
in  Chapter  V. 
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II.  Backround 


Although  the  breakdown  of  gases  has  been  studied  ex- 
tensively for  many  years,  gas  breakdown  under  the  action 
of  radiation  at  optical  frequencies  is  a more  recently  ob- 
served(1962)  phenomenon.  Since  fluxes  intense  enough  to 
produce  breakdown  at  optical  frequencies  can  only  be  pro- 
duced by  a laser,  it  is  within  this  context  that  the  phe- 
nomenon is  studied.  The  theory  involved  has  been  firmly 
established  and  will  be  discussed  in  this  chapter.  The 
first  part  of  this  chapter  is  devoted  to  a qualitative 
discussion  of  the  processes  involved,  while  the  remainder 
of  the  chapter  will  deal  with  some  of  the  major  theories 
utilized  to  explain  breakdown. 

Breakdown  (A  Qualitative  Discussion) 

The  breakdown  of  a gas  will  occur  if,  during  the  time 
of  the  laser  pulse,  the  gas  can  absorb  enough  energy  to 
achieve  and  maintain  a high  degree  of  ionization.  Pro- 
cesses which  could  account  for  this  ionization  are  the  ab- 
sorption of  photons  by  a neutral  atom,  inelastic  colli- 
sions between  atoms  and  electrons,  and  a combination  of 
both  processes.  Those  processes  which  could  reduce  the 
degree  of  ionization  are  the  recombination  of  electrons 
and  ions,  and  the  diffusion  of  the  charged  particles  from 
the  volume  of  the  gas  upon  which  the  laser  radiation  is 
incident. 
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The  direct  ionization  of  atoms  by  laser  photons  was 
studied  extensively  by  Keldysh(Ref  1).  Since,  in  most 
cases,  the  value  of  the  ionization  energy (5-20  eV)  is  much 
greater  than  the  photon  energy (0.1-1. 5 eV)  it  is  impossi- 
ble for  one  photon  to  ionize  the  atom.  However,  if  I/hw 
photons  can  be  absorbed  simultaneously  by  an  atomic  elect- 
ron, ionization  would  be  possible.  At  atmospheric  pres- 
sure the  laser  intensities  required  to  produce  breakdown 
by  this  mechanism  are  much  greater  than  those  under  which 
breakdown  is  observed  experimentally.  Thus,  the  "multi- 
quantum photoeffect"  is  not  considered  to  be  a likely 
breakdown  mechanism. 

The  process  of  electron  cascade  ionization(Ref  2) 
has  been  accepted  as  the  mechanism  most  likely  involved 
in  breakdown.  The  cascade  can  begin  with  a single  "prim- 
ing electron".  This  electron  absorbs  light  quanta  during 
collisions  with  atoms.  Along  with  this  inverse  bremsstrah- 
lung  absorption  of  energy  by  the  electrons  there  occurs  a 
stimulated  bremsstrahlung  emission.  The  absorption,  how- 
ever, is  more  likely  than  emission  so  that  over  a period 
of  time  a net  absorption  of  energy  occurs.  This  net  ab- 
sorption of  energy  can  be  thought  of  as  a diffusion  of  the 
electron  along  a one-dimensional  energy  axis.  If  excita- 
tional  processes  are  temporarily  ignored,  the  energy  of  tne 
electron  will  continue  to  increase  until  it  exceeds  the 
ionization  potential  of  the  neutral  atoms.  It  will  then 
collide  inelastically  with  one  of  the  atoms  and  ionize  it. 
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The  result  will  be  two  electrons  of  very  low  energy  re- 
placing the  one  electron  with  energy  slightly  greater  than 
I.  These  two  electrons  can  now  absorb  energy  and  repeat 
the  process.  This  cascade  will  continue  until  breakdown 
is  achieved. 

Energetic  electrons,  however,  are  also  capable  of 
causing  excitations.  Since  excitation  requires  less  en- 
ergy than  ionization  an  electron  may  lose  its  energy  to 
excitation  many  times  before  acquiring  an  energy  suffi- 
cient for  an  ionization.  Under  conditions  of  sufficiently 
high  flux,  atoms  excited  into  higher  atomic  states  can  be 
rapidly  ionized  by  the  absorption  of  two  or  three  quanta. 
In  this  case  an  excitation  is  equivalent  to  an  ionization 
and  the  cascade  continues  unhindered.  At  lower  fluxes  and 
photon  energies  of  less  than  about  1 eV,  photoionization 
of  the  excited  states  becomes  unlikely(for  the  same  reason 
as  multiphoton  ionization)  and  energy  going  into  excita- 
tions slows  the  development  of  the  cascade.  Once  an  atom 
is  excited  the  probability  that  it  will  encounter  a second 
electron  and  be  ionized  is  extremely  small  except  during 
the  latter  stages  of  the  cascade  when  the  density  of  elec- 
trons and  excited  atoms  is  high(Ref  3:658). 

The  development  of  the  cascade  can  also  be  hindered 
by  diffusion  of  the  free  electrons  from  the  region  of  ac- 
tion of  the  incident  radiation.  The  diffusion  losses  will 
be  smaller  at  higher  pressures  and  larger  beam  diameters. 
At  fluxes  of  interest(S  > lO^W/cm2),  diffusion  losses  can 
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be  assumed  to  be  negligible  at  pressures  greater  than  1 
atm  and  beam  diameters  larger  than  10  cm(Ref  2:773). 
Another  mechanism  that  will  hinder  the  cascade  is  three 
body  recombinations.  This  mechanism,  however,  is  not  con- 
sidered significant  since  losses  due  to  it  will  only  be 
observed  very  late  in  the  cascade. 

Since  at  least  one  electron  must  be  present  for  the 
cascade  to  proceed,  the  question  of  the  origin  of  the  first 
electron  is  often  raised.  The  presence  of  a small  quantity 
of  easily  ionizable  impurity  atoms  is  one  source.  The  pro- 
duction of  the  first  electron  can  then  occur  as  a result 
of  photoionization  or  collision  involving  the  impurity 
atom.  Another  source  of  electrons  is  the  atmosphere.  Non- 
equilibrium atmospheric  effects  such  as  cosmic  radiation 
result  in  an  average  free  electron  density  of  100/cm  . 
Therefore,  under  appropriate  conditions,  the  first  elec- 
tron can  be  found. 

Finally,  if  R is  defined  to  be  the  cascade  growth 
rate  then  the  number  of  free  electrons  at  any  time  t can 
be  expressed  by 


n 


e 


(1) 


where  n°  is  the  initial  free  electron  density.  Now,  if 
e 

the  breakdown  criterion  is  chosen  to  be  the  achievement 
of  a specific  electron  population  ne,  then  the  time  re- 
quired for  breakdown  to  occur  is  given  by 
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(2) 


*1  =Tln(ne/ne) 


where  r=  1/R  is  the  cascade  time  constant.  As  can  be 
seen  from  Eq  (2)  the  occurrence  of  breakdown  is  not  strong- 
ly dependent  on  the  initial  electron  density.  It  is,  how- 
ever, highly  dependent  on  the  time  constant  r which  is  in 
turn  a strong  function  of  the  laser  flux.  It  is  for  this 
reason  that  an  abrupt  breakdown  threshold  is  found  to 
exist(Ref  5:657)* 

The  Boltzmann  Equation 

Early  investigations  revealed  that  the  problem  of 
laser  induced  gas  breakdown  could  be  best  approached  through 
the  use  of  the  Boltzmann  equation 

II  + V-Vt  + I vTf  ■=  g£|ooll  (5) 

where  f is  the  electron  distribution  function,  v the  elec- 
tron velocity, 

represents  the  change  in  f due  to  collisional  processes. 
Under  the  assumption  of  spatial  homogeneity  of  f the  se- 
cond term  in  Eq  (3)  vanishes. 

In  a manner  identical  to  that  used  for  microwave 
breakdown(Ref  4),  an  electron  heating  term  can  be  derived. 
This  term  describes  the  rate  at  which  an  electron,  acted 
upon  by  a field,  gains  energy  on  collision  with  an  atom. 


a the  electron  acceleration,  and 


*£| 


coll 
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Under  these  conditions  Newton's  equation  of  motion  for  the 
electrons  can  be  written  as(Ref  4:170) 


m ~ + m v v = e E e-^**^ 
e dt  e 


(4) 


where  v is  the  electron-atom  collision  frequency,  e the 
electronic  charge,  and  w is  the  angular  frequency  of  the 
applied  field,  E.  Solving  Eq  (4)  for  v yields 


v = 


eEe~iwt 

m v -iwj 
0 


(5) 


The  third  term  of  Eq  (3)  can  be  rewritten  as 


a 


(mev) 


(6) 


so  that  with  a = ~ e_1  and  v given  by  Eq  (5)»  Eq  (6) 

me 

averaged  over  a period  becomes 


S-vv* 


e2E2v  5f 


(7) 


This  equation  can  also  be  derived  in  a more  rigorous  man- 
ner directly  from  the  Boltzmann  equation (Ref  4:46-47,166). 
Substituting  Eq  (7)  into  Eq  (3)  yields  a final  form  for 
the  Boltzmann  equation 


n_  dt 


e 


^Icoll 


(8) 
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where  f is  defined  to  be  a function  of  energy  and  i is 
given  by 


e 


n e S v 
e 

* — *— 

mec  «0(v  +w  ^ 


(9) 


p 

where  the  substitution  S = («q/2)E  c was  made  in  Eq  (7),  S 

2 

being  the  laser  flux  in  W/cm  . Included  in  the  collision 
term  on  the  right  side  of  Eq  (8)  are  the  rates  of  excita- 
tion and  ionization  as  well  as  the  rates  of  any  other  sig- 
nificant collisional  processes. 


A Temperature  Model 

If  the  collisional  terms  are  written  out  explicitly 
Eq  (8)  becomes 


§LL±1  + _L  2£lL 1 

Ai.  ' ~ _ 


at 


ne  dt 


\ ?>ini^(<)f(°  + ± X niR^(€+  xi;j)f(€  + x.d) 

- £ n .R*(«)f  (O  + £ ndR*( t + Ij)f ( « + Ij) 

<3  • «J 

♦ 1(0  I nj  /RJ(<)f(«)d<  + R (10) 

j 

where  R^j  and  Rj  are  the  rates  for  excitation  and  ioniza- 
tion, and  Ij  represent  the  excitation  and  ionization 

energies,  and  the  ni  are  the  population  densities  of  the 
various  atomic  levels.  The  first  four  terms  on  the  right 
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side  of  Eq  (10)  represent  the  rate  at  which  electrons  of 
energy  c are  either  gained  or  lost  to  the  distribution  as 
a result  of  excitation  and  ionization  of  atoms.  The  fifth 
term  represents  the  rate  at  which  electrons  of  zero  energy 
are  created  by  ionization.  The  last  term,  I?,  represents 
all  reverse  rates(deexcitations,  recombinations)  and  is 
assumed  to  be  small  during  the  early  stages  of  breakdown 
if  most  of  the  atoms  are  originally  assumed  to  be  in  a 
ground  state. 

Now,  if  the  gas  is  assumed  to  be  made  up  of  atoms 
which  possess  a nucleus  and  two  electronic  levels,  a ground 
state  and  an  excited  state,  Eq  (10)  simplifies  to 

+ SLLll  m 

ne  d( 

-n1Rx(«)f(<)  + n1Rx( < + x)f( « + x) 

-n1RI(«)f(<)  + n1RI( € + I)f(  « + I) 

+ *(<)/niRI(<)  f(«)  d«  (11) 

At  this  point,  the  exact  form  of  the  electron  distribution 
function  is  unknown.  If  it  is  assumed  to  be  maxwellian, 

«•>  - 2 #§  °-‘/T  <12) 
where  T is  the  electron  temperature,  then  a temperature 
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model  approach  towards  breakdown  will  emerge.  Taking  the 
zeroth  energy  moment  of  Eq  (11)  yields  an  equation  for  elec- 
tron growth (Appendix  C) 

dn  T 

3F~  = <R  > nen1 

where  <PX>,  the  ionization  rate,  is  a function  of  temper- 
ature, obtained  by  averaging  R1  over  all  energy.  Taking 
the  first  energy  moment  yields  an  energy  balance  equation 
(Appendix  C) 

^ - £-  - *<«*>  ”1  - «rI>  "i 

e 

dn 

- m—Z  (14) 

ne  dt 

where  <c>  , the  average  energy  in  eV,  is  given  by  ^ T. 

Eqs  (13)-(14)  taken  together  can  be  used  to  investigate 
and  predict  many  of  the  phenomena  involved  in  breakdown 
and  are  often  referred  to  as  a "temperature  model". 

The  conditions  required  for  breakdown  csui  now  be  es- 
timated using  Eqs  (13)-(14).  Under  the  assumption  that 
the  average  electron  energy  equilibrates  rapidly (t  1**) 

Eq  (14)  becomes 

- x<R*>  + I<RX>  n^  + <€XRI>  (15) 


where  Eq  (13)  was  also  used.  The  rates  can  be  expressed 


as (Ref  5:388-390) 


<RX>  . Rx  e-x/T 

(16) 

and 

<BJ>  = hJ  e*I/T 

(17) 

y T m 

where  Rq  and  Rq  are  slowly  varying  functions  of  T which 

x X 

can  be  assumed  constant.  For  R,»  R^  the  last  two  terms  on 

o o 

the  right  side  of  Eq  (15)  can  be  ignored.  Solving  Eq  (15) 
for  e_x yields 

e-x/T  , (18) 

nen1x  Ro 

Now,  if  both  sides  of  Eq  (18)  are  raised  to  the  I/x  power 
and  the  results  substituted  into  Eq  (17) » <R*>  becomes 


(19) 


Solving  Eq  (13)  for  n and  applying  the  breakdown  criterion 

© 

yields 


ln(n2/n°) 


n1Bo 


e'  e‘ 


n 


e2S  v 


1X  *omec  o(^ 


5 


I/x 


(20) 


which  is  Eq  (2)  with  r = 1/n/)<R^>,  and  replaced  by  Eq  (9). 
Thus,  for  a given  laser  intensity  S,  breakdown  will  occur 
only  if  the  laser  pulse  length  meets  the  criterion  of 
Eq  (20). 
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The  condition3  under  which  the  assumption  of  a temp- 
erature is  valid  will  be  discussed  later  in  this  report 
(Chapter  IV.) 

The  Quantum  Kinetic  Equation 

The  classical  electron  heating  rate  of  Eq  (9)  was 
originally  derived  for  use  in  microwave  breakdown  studies 
sind  assumes  a continuous  absorption  of  energy  by  the  elec- 
tron. At  microwave  frequencies  this  assumption  is  valid 
since  the  energy  of  each  photon  is  small,  and  the  average 
energy  per  collision  absorbed  by  an  electron(as  calculated 
by  Eq  9)  is  many  times  the  energy  of  a single  photon.  At 
optical  frequencies  the  energy  of  a photon  becomes  appre- 
ciable and  the  average  energy  per  collision  absorbed  by  an 
electron  is  only  a small  fraction  of  the  energy  of  a sin- 
gle photon.  Thus,  an  electron  will  experience  many  colli- 
sions which  do  not  involve  the  absorption  of  energy  and  a 
quantum  mechanical  approach  should  be  taken(Ref  3:657)* 

The  heating  of  electrons  by  the  radiation  is  the  net 
result  of  two  competing  processes.  The  first,  bremsstrah- 
lung  absorption,  occurs  when  an  electron  absorbs  a photon 
from  the  radiation  field  during  collisions  with  atoms  or 
ions.  The  second,  bremsstrahlung  emission,  occurs  when  an 
electron  is  stimulated  by  the  radiation  field  to  emit  a 
photon  upon  collision  with  an  atom  or  ion.  A coefficient 
of  bremsstrahlung  absorption  can  be  derived  quantum  mecha- 
nically and  is  given  by(Ref  6:540) 
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a(.)  - 2e2  s\2^*  *%<**  W2)  (21) 

3mec<ow  «-  e J 

where  « is  the  electron  energy,  ftw  is  the  energy  of  the 
photon  absorbed,  and  om  is  the  momentum  transfer  cross 
section  of  the  electron-atom  pair.  A coefficient  of  brem- 
sstrahlung  emission  cam  be  obtained  from  Eq  (21)  using  mic- 
roreversibility (Ref  2:774) 


b(  * + ft-)  = a(«) 


r_i_T 

(«  + ft«J 


(22) 


Now,  if  the  quantum  mechanical  rates  for  bremsstrahlung 
emission  and  absorption  are  used  to  replace  the  microwave 
heating  term  in  Eq  (10)  the  resulting  quantum  kinetic 
equation  will  be 


a|a(«-  ft-)f(«-  ft-)  - (a(«)  + b(«)]f(«) 
+ b(«+  ft-)f(«  + ft-))  - X niR^(«)f(«) 

+ i,f>iniR^(‘+  ’‘U5  * £ ndHjCOf(,) 

+ I n;jRj(«+  Ij)f(«+  Ij) 

J 

+ 4(0  X n 7R*(0f(0<i« 

3 d*  0 


(23) 


where  F is  the  photon  flux(#/cm-sec)  and  n is  the  total 
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number  of  atoms.  The  first  term  on  the  right  hand  side  of 
Eq  (23)  is  the  rate  at  which  electrons  of  energy  < are 
gained  as  a result  of  bremsstrahlung  absorption  and  emis- 
sion. 

A Diffusion  Model 

As  previously  mentioned  the  heating  of  the  electrons 
by  a laser  can  be  viewed  as  a diffusion  of  electrons  along 
the  energy  axis.  This  approach  was  originally  presented 
by  Zeldovich  and  Raizer(Ref  2).  Eq  (23)  can  be  rewritten 
in  the  form 

* Fnfl{a(«-  hw)f(«-  h«)  - |a(0  + b(«)]f(«) 

+ b(«+  fi**)f(«+  tw)}  + Q (24) 

where  Q includes  all  electron  energy  gain  and  loss  mecha- 
nisms not  specifically  written  out.  In  the  region  near 
the  excitation  and  ionization  energies  hw/«  is  a small 
quantity  and  the  bremsstrahlung  absorption  and  emission 
terms  cam  be  expanded  in  a Taylor  series  in  powers  of 
h»/«  . Retaining  terms  to  second  order  results  in 
(Ref  2:775) 

♦ Q (25) 
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J = - D + uf(«) 

(26) 

where 

D = ■£  Fna(hw)2[a(0  + b(«)] 

(27) 

and 

u = Fna&a[a(0  - b(«)]  - ^ 

(28) 

Eqs  (25)-(26)  represent  a diffusion  of  electrons  along  the 
energy  axis.  D is  the  diffusion  coefficient,  u is  the  net 
velocity  of  the  electrons  along  the  energy  axis,  and  Q re- 
presents an  electron  "source".  Eqs  (27)-(28)  can  be  sim- 
plified further  by  expanding  b(«)  in  terms  of  hw/« 
yielding(Ref  2:775) 


D = Fn_(hw)2a(«) 

51 


(29) 


(50) 


If  Q is  assumed  to  be  zero  (no  sources),  Eq  (25)  can  be 
multiplied  by  « and  integrated  over  all  energy  so  that  in 
the  limit  as  Iiu  /<-»  0 it  becomes(Ref  2:776) 


d <*> 

“3t“ 


e2S  v 


me  < w 
e o 


(51) 


where  v was  defined  to  be  n_vo ■ . At  optical  frequencies 

CL  ID 

v2«  «2,  and  Eq  (31)  agrees  with  Eq  (9)  which  gives  the 
classically  defined  electron  heating  rate.  Thus,  it  can 
be  seen  that  the  quantum  kinetic  model  reduces  to  the 
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correct  classical  limit  as  hwA-^O,  and  for  Hw/t  small 
the  classical  heating  rate,  t , can  be  used  to  accurately 
describe  the  net  bremsstrahlung  absorption  of  laser  radia- 
tion. 

The  exact  method  of  solution  of  Eq  (25)  is  not  impor- 
tant, but  the  results  does  have  a natural  physical  inter- 
pertation  which  should  be  discussed.  Q,  in  Eq  (25),  can  be 
replaced  by  a boundary  condition  involving  the  probability 
of  ionization.  The  electrons  are  then  assumed  to  increase 
in  energy  until  they  reach  IQ,  an  energy  1-3  eV  greater 

a 

than  the  ionization  potential.  Having  reached  this  energy 
the  electrons  will  lose  their  energy  to  ionization  with  a 
probability  a and  to  excitation  with  a probility  1-  a . 

Under  these  conditions  Eq  (25)  can  be  solved  for  the  cas- 
cade time  constant (Ref  2:777) 


t = t /a  (32) 

cl 

where  t = I /3u  is  the  time  it  takes  an  electron  to  grow 
a a 

in  energy  from  0 to  I , and  1/«  represents  the  number  of 
times  this  cycle  must  be  repeated  before  sin  ionization 
occurs.  Thus,  the  electron  growth  rate,  R = 1/t,  depends 
on  the  laser  flux  through  t and  on  the  properties  of  the 
gas  through  «.  This  is  consistent  with  the  results  of  the 
temperature  model. 
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III.  Approach 


If  a comparison  of  the  results  of  a temperature  model 
with  those  of  the  more  accurate  quantum  kinetic  model  is 
to  be  valid  care  must  be  taken  to  insure  that  each  model  is 
solved  under  similar  conditions.  The  first  part  of  this 
chapter  discusses  the  various  assumptions  made  to  develop 
self  consistent  models.  The  remainder  of  the  chapter  out- 
lines the  solution  of  the  models  and  the  comparisons  made. 

The  Model  Atom 

The  gas  upon  which  the  laser  energy  is  incident  con- 
sists of  a single  atomic  species.  The  atom  has  a mass  of 
15  amu  and  consists  of  a core  and  one  electron.  Three 
energy  levels  are  permitted:  a ground  state  at  -14-  eV,  a 
first  excited  state  at  -3.5  eV,  and  a second  excited  state 
at  -I.54-  eV.  Ionizations  can  occur  from  any  state  and 
transitions  are  allowed  between  any  two  levels.  Multiple 
ionizations  of  an  atom  are  ignored.  The  temperature  of 
the  gas  (atoms  and  ions)  is  assumed  to  remain  constant 
until  late  in  the  cascade  when  heating  by  the  electrons 
becomes  significant. 

Cross  Sections  and  Rates 

Cross  sections  and  rates  for  the  various  processes 
involved  in  breakdown  were  either  chosen  or  calculated  so 
as  to  insure  complete  correspondence  between  the  quantum 
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kinetic  rate  and  its  temperature  model  counterpart . The 
quantum  kinetic  rates  for  ionization  and  excitation  were 
determined  using  analytic  formulas 


R(‘)  *(0 


(33) 


where  <r(«)  is  the  cross  section  for  either  ionization  or 
excitation.  Substituting  the  ionization  cross  section 
into  Eq  (33)  yields (Ref  5:392) 


hJo  . 


j£-g  1 I,) 


«2I. 

J 


(34) 


where  is  the  ionization  rate  for  the  jth  level  and  I. 

o J 

is  the  corresponding  ionization  energy.  In  a similar  man- 
ner an  excitation  rate  can  be  obtained(Ref  5:396) 


Rid(0  ■ 777  *i;i? V3/2(<-  V *iXio  (35) 


where  R*^  is  the  excitation  rate  from  level  i to  level 

x . . is  the  excitation  energy,  and  f.  . the  oscillator 
ij  ij 

strength  for  absorption.  As  can  be  seen  from  Pig  1 these 
analytic  rates  compare  favorably  with  similar  data  obtained 
experimentally(Ref  7:'I15)« 

The  corresponding  temperature  model  rates  were  obtain- 
ed by  multiplying  the  quantum  kinetic  rates  by  a Maxwell- 
Boltzmann  distribution  function  and  integrating  over  all 
energy.  The  ionization  rate  is 
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rJ(T) 


,3/2^  4 


VBT(4 


■j£T3 T7?  [e*I3/T  + Ei('1.1/T)]  <36) 

w«  ) I . T ' 1 J 

A ' 1 * 


e'  o'  o 


where  is  a tabulated  integral  defined  by 


VJt 

E^(x)  * J dt 

-oo 


(37) 


Similarly , the  excitation  rate  is  given  by 

P3/2_  4, 

RX  (T)  - 2 **  2 7n  re”xii/T 

RiJU;  L J 

^(4w<o)  xij  T 


+ ' ^ Ei(-xid/T)] 


(38) 


A complete  treatment  of  the  above  cross  sections  and  rates, 
along  with  a table  of  numerical  formulas,  can  be  found  in 
Appendicies  A and  B.  Rates  for  deexcitation  and  recombi- 
nation were  not  computed  since  they  were  ignored  in  the 
solution  of  the  equations. 

Finally,  a value  for  the  electron-atom  momentum  trans- 

o 

fer  cross  section,  a-  = 7«1  x 10  cm  , was  taken  from 

* m 

Brown(Ref  7:5)»  This  cross  section  is  used  to  compute  the 
value  of  v in  Eq  (9)  and  also  to  compute  a(«)  in  Eq  (21). 

A cross  section  for  electron-ion  collisions  was  not  cal- 
culated since  electron-ion  collisions  were  ignored. 
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The  Quantum  Kinetic  Model 

For  the  three  level  model  atom  assumed  the  quantum 
kinetic  equation  can  be  written  as 


■ Fna{a(«-  hw)f(«-  h«)  - [a(«)  + b(«)]f(«) 


+ b(«  + hw)f(«+  hu.)}  - r n.R^.(Of(«) 

i,  j>i  1 


+ i,iiniR^(<+  Xid)f(<+  Xij}  ■ J ndV°f(0 


+ L njRj(<+  Ij)*(«  + 


+ 1(0  L n.  /Rj(* )f  (« )d< 

-5  JO  J 


(23) 


The  first  term  on  the  right  side  of  Eq  (23)  represents  the 
rate  at  which  electrons  of  energy  « are  gained  as  a result 
of  bremsstrahlung  absorption  and  emission.  The  coeffi- 
cients a and  b are  given  by  Eqs  (21)-(22);  F is  the  photon 
flux(#/cm  -sec)  and  n is  the  total  number  of  atoms 

gL 

(#/cm^).  The  next  two  terms  of  Eq  (23) , involving  R*^, 
represent  the  rate  at  which  electrons  of  energy  t are 
lost  as  a result  of  excitations.  The  summation  over  i 
goes  from  1 to  3 and  the  R*^  are  given  by  Eq  (35).  The 
fifth  and  sixth  terms  of  Eq  (23)  represent  the  rate  at 
which  electrons  of  energy  < are  lost  as  a result  of 
ionizations.  The  summation  over  j goes  from  1 to  3 and 
the  Rj  are  given  by  Eq  (34) • Finally,  the  last  term  of 
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Eq  (23)  represents  the  rate  at  which  electrons  with  zero 
energy  are  created  by  ionizations  from  each  of  the  three 
atomic  levels. 

There  are  several  assumptions  implicit  in  Eq  (23) 
that  should  be  mentioned.  First,  the  diameter  of  the 
incident  laser  beam  is  assumed  to  be  large(~1  cm). 
Diffusion  of  electrons  from  the  region  of  the  beam's 
action  will  then  be  minimal ( Tj~10“^sec)  and  can  be 
neglected.  During  the  early  portion  of  the  cascade, 
electron-ion  collisions  will  be  rare  and  can  be  neg- 
lected. This  eliminates  the  need  to  consider  the 
bremsstrahlung  heating  of  electrons  by  ions,  and  triple 
body  recombination  resulting  from  the  collision  of  an 
ion  and  two  electrons.  For  a similar  reason,  electron- 
electron  collisions  can  be  ignored.  Finally,  the  deexci- 
tation of  atoms  is  neglected.  The  energy  lost  to  deexci- 
tations is  very  small  and  the  rate  of  ionization  from  the 
excited  states  is  not  strongly  affected. 

Besides  Eq  (23),  three  other  equations  contribute 
to  the  quantum  kinetic  description  of  breakdown.  These 
are  the  equations  governing  the  change  in  the  popula- 
tion of  the  three  atomic  levels 

- -n1ne[<R^(«)>  + <R*2(«)>  + <R^(0>]  (39) 

||2  - n1ne<R^f2(«)>  - n2ne[<R^(0>  + <R!3(«)>]  (40) 


24 


||3  = nine<R*3(0>  + n2ne<R*3(0>  - n3ne<R*(«)> 


(41) 


where  <R(«)>  is  the  energy  average  of  the  corresponding 
rate.  These  equations  together  with  Eq  (23)  can  now  be 
solved  numerically  using  methods  similar  to  those  of  Ref  8 
(see  Appendix  D).  It  should  be  noted  that  a good  estimate 
of  the  cascade  growth  rate  can  be  obtained  by  solving  Eq 
(23)  independent  of  Eqs  (39)-(41)  using  the  initial  values 
of  n^,  n2,  n3«  However,  Eqs  (39)-(41),  as  well  as  other 
effects,  become  important  if  a solution  is  to  be  carried 
to  breakdown. 

The  Temperature  Model 

In  an  effort  to  insure  self-consistency,  the  temp- 
erature model  equations  are  obtained  directly  from  the 
quantum  kinetic  equations.  For  ftuA  assumed  small,  the 
bremsstrahlung  absorption  and  emission  terms  in  Eq  (23) 
can  be  replaced  by  the  microwave  heating  term,  Eq  (9)* 

The  electron  distribution  function,  an  unknown  in  Eq  (23) » 
is  assumed  to  be  maxwellian 

e“/T  (42> 

If  Eq  (23),  as  modified,  is  integrated  over  all  energy 
the  zeroth  energy  moment  of  the  Boltzmann  equation  will 
result 
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where  the<Rj>are  now  averages  over  a maxwellian,  and 
are  a function  of  temperature.  Eq  (43)  describes  the  rate 
of  growth  of  free  electrons.  Eq  (23)  can  also  be  multi- 
plied by  < and  integrated  over  all  energy  to  obtain  the 
first  energy  moment  of  the  Boltzmann  equation 


d«> 

“3F“ 
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L n • x . .(R.  - 1 -vt— 

i,d>i  1 10  10  ne  ^ 


(44) 


where  «>=  (3/2)T  is  the  average  electron  energy.  Eq  (44) 
describes  the  rate  of  growth  of  the  average  free  electron 
energy.  To  complete  the  description  of  the  temperature 
model  Eqs  (39)-(41)  must  be  included  with  Eqs  (43)-(44). 
The  rates  in  these  equations  are  all  functions  of  temper- 
ature and  are  the  result  of  averaging  over  a maxwellian. 

A more  complete  derivation  of  the  temperature  model  can 
be  found  in  appendix  C. 

The  assumptions  involved  in  the  temperature  model  are 
identical  with  those  of  the  quantum  kinetic  model  with  two 
exceptions.  The  bremsstrahlung  emission  and  absorption 

terms  were  replaced  by  a microwave  heating  term.  For  fiw/* 

* 

small  the  two  terms  will  yield  the  same  electron  heating 
rate  so  that  consistency  is  maintained.  The  second  excep- 
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tion  involves  electron-electron  collisions.  If  the  max- 
wellian  nature  of  the  electron  distribution  is  to  be  main- 
tained, electron-electron  collisions  must  be  at  least  im- 
plicitly allowed.  Since  it  is  the  validity  of  the  assump- 
tion of  a maxwellian  that  is  being  tested,  this  will  not 
cause  a conflict. 

The  temperature  model  equations  are  solved  numeri- 
cally using  a procedure  outlined  in  appendix  E. 

Solution  and  Comparison 

The  two  breakdown  models  have  been  solved  under  a 
variety  of  conditions  to  obtain  a comparison  of  the  re- 
sults and  behavior  of  each  model.  Certain  initial  condi- 
tions were  imposed  and  remained  fixed  for  all  comparisons. 
The  laser  radiation  was  assumed  to  have  a wavelength  of 
10.6K  . This  fixes  the  energy  of  a photon  at  .12  eV.  The 
initial  density  of  ground  state  atoms  was  assumed  to  be 
2.5*  10  Vcm  . The  initial  temperature  of  the  gas  was  set 
at  .1  eV.  This  condition  fixed  the  initial  population  den- 
sity of  the  first  and  second  excited  states  as  well  as  the 
total  number  of  atoms  initially  present.  The  cascade  was 
assumed  to  proceed  from  an  initial  free  electron  density  of 
100/cm^.  Since  solution  of  the  quantum  kinetic  model  re- 
quires that  the  energy  axis  be  divided  into  bins (see  ap- 
pendix D),  a standard  of  25-1 eV  bins  was  chosen. 

Under  these  conditions  the  two  models  were  solved 
utilizing  the  rates  found  in  table  X(Appendix  B).  The 
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cascade  growth  rates  predicted  by  each  model  were  computed 

7 'I  P 

for  a wide  range  of  incident  flux  densiti es(5  x ICr  - 10 

p 

W/cm  ) and  compared.  A study  of  the  apportionment  of  en- 
ergy among  the  various  loss  mechanisms  suggested  that  a 
comparison  of  the  growth  rates  predicted  by  the  two  models 
be  performed  assuming  a two  level  atom.  Finally,  the  re- 
sults of  the  quantum  kinetic  model  were  compared  for  vari- 
ous bin  widths. 

To  determine  the  effects  of  manipulation  of  energy 
levels  and  ionization  and  excitation  rates,  comparisons  of 
the  growth  rates  were  made  under  the  following  conditions: 

1.  increased  by  a factor  of  10. 

2.  decreased  by  .5 

3.  R^"  increased  by  a factor  of  10. 

4.  R*  decreased  by  .5 

5.  Ground  state  energy  raised  to  -13  eV. 

6.  Ground  state  energy  lowered  to  -15  eV. 

7.  1st  excited  state  energy  raised  to  -2.5  eV. 

8.  1st  excited  state  energy  lowered  to  -4.5  eV. 

The  results  of  all  comparisons  performed  will  be 

discussed  in  Chapter  IV. 

An  Improved  Diffusion  Model 

In  an  attempt  to  find  a breakdown  model  applicable 
over  a wider  range  of  incident  fluxes,  the  improved  dif- 
fusion model  of  Vyskrebentsev  and  Raizer(Ref  9)  was  stu- 
died. Unlike  the  old  model  of  Zeldovich  and  Raizer, 
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excitations  are  treated  specifically  rather  than  being 
lumped  into  a generalized  term,  Q.  Assuming  a two  level 
atom,  Eq  (25)  can  be  modified  to(Ref  9:55) 

- - lr  - b(«)f(«)v*  (45) 

(0  for  «<x 

a (46) 

1 for  x<  «<  I 
a a 

The  energy  x_  is  defined  to  be  approximately  leV  greater 
than  the  excitation  energy,  x,  while  I is  leV  greater 
than  the  ionization  energy,  I.  In  the  interval  from  x„ 

6L 

to  I_,  the  average  excitation  rate,  v*,  is  assumed  con- 

cl 

stant. 

Eq  (45)  can  be  solved  with  the  aid  of  several  bound- 
ary conditions.  At  « = x_,  f(0  and  df(«)/d<  are  assumed 

OL 

to  be  continuous.  All  electrons  reaching  I are  assumed  to 
lose  their  energy  to  ionization  with  a probability  0 and 
to  excitation  with  a probability  (1-0)  so  that  f(IQ)  = 0. 
Finally,  the  rate  at  which  electrons  appear  with  energy 
zero  equals  the  rate  of  excitation  in  the  interval  x<«<I, 
plus  twice  the  rate  of  ionization(this  accovints  for  newly 
born  electrons).  Under  these  conditions  Eq  (45)  yields  a 
transcendental  equation  for  the  cascade  growth  rate,  v^(Ref 
9:54), 
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e^-^tcosMf)  + 2 sinh(f)]-  e'^-1  [coah(|) 


- z sinh(^)]  - 2a(1+0)y  - 2(1-z“2){  y cosh  [(a-1)y] 


where 


+ sinh[(a-1)y]  - ay}  = 0 


y = vfe(  v±  + v*)xa/(e/ne) 


z = VI  + v*/v ^ 


a = Va 


(47) 

(48) 

(49) 

(50) 


Before  Eq  (47)  can  be  solved  0,  v*,  and  e must  be  evaluated. 
The  value  of  P is  obtained  from 


P 


R^«a)  ♦ R^da) 


(5D 


and  is  found  to  be  .09  . The  rates  in  Eq  (51)  were  calcu- 
lated using  table  X(Appendix  r-).  The  value  v*  is  defined 
by (Ref  9:53) 


V*  ■ nav*** 


(52) 


where  v*  and  7*  are  mean  values  of  the  electron  velocity 
and  excitation  cross  section  in  the  interval  xaf:<<Ia 
determined  by 
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for  n = 2.5  x 10^/cm^  Eqs  (52)-(54)  yield  a value  for  v* 

cl 

AA 

of  8.5  x 10  /sec.  The  quantity  I:  is  evaluated  over  a 
range  of  fluxes  using  Eq  (9)*  Eq  (47)  can  then  be  solved 
numerically  and  the  growth  rates  compared  with  correspon- 
ding quantum  kinetic  rates. 


A Two-Temperature  Model 

Analysis  of  the  results  of  the  comparisons  between 
the  temperature  model  and  the  quantum  kinetic  model (Chap- 
ter IV)  suggested  that  modification  of  the  temperature 
model  distribution function  might  lead  to  a more  accurate 
prediction  of  the  growth  rate.  The  most  encouraging  ap- 
proach appeared  to  be  one  in  which  the  electron  distribu- 
tion function  was  broken  into  two  sections,  each  described 
by  a different  temperature. 

The  problem  is  first  simplified  by  considering  a two 
level  atom.  The  two-temperature  electron  distribution 
function  is  then  defined  as 
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for  0<€<x 


/C,  ‘/T1 

f(0=<  (55) 

lc2</1>/2e“  <//T2  for  xsk- 

where  x is  the  excitation  energy  of  the  model  atom.  This 
distribution  function  is  assumed  to  be  continuous  at  the 
boundary  of  the  two  regions 

C1e_x/T1  = C2e“x/T2  (56) 

A second  condition  requires  that  the  distribution  function 
be  normalized 

0.  /«1/2e"  ‘/T1  d€  + Cp  7«1/2e~*/T2  d«  « 1 (57) 

1 0 d x 

The  energy  balance  equation  of  the  old  temperature  model 
is  also  needed 

- *<r1>”i  - sf 

e e 

It  should  be  noted  that  the  average  energy  and  the  rates 
in  this  equation  must  be  reevaluated  using  the  modified 
distribution  function  defined  in  Eq  (55)-  Assuming  a 
rapid  equilibration  of  the  average  electron  energy  and 
substituting  for  dn  /dt  using  Eq  (13)  the  energy  balance 
equation  becomes 
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(58) 


— - xn1<Rx>  - (I  + <«>)n1<RI>  = 0 
e 

Eqs  (56)-(58)  form  a set  of  three  equations  in  four 
unknowns.  A fourth  condition  on  the  unknowns  can  be  ob- 
tained by  taking  the  second  energy  moment  of  the  Boltzmann 

equation (Appendix  F).  This  will  yield  an  equation  for  the 

2 

time  rate  of  change  of  « > 

= 2^<«>  + x2n/J<Rx>  + I2n/J<RI> 
e 

_ 2 dn 

- 2xn1<«Rx>  - 2In1«R1>  = ^ (59) 

e 

where  the  averages  must  be  taken  over  the  distribution 
function  defined  by  Eq  (55) • Again,  assuming  a rapid 

p 

equilibration  of  <«  > and  using  Eq  (13)  yields 

2-^-<«>  + x2n1<Rx>  + (I2  - <*2>)n1<RI> 
ne 

-2xn/j<*Rx>  - 2In^<«R^>  = 0 (60) 

The  two-temperature  model  consisting  of  Eqs  (56)-(58) 
and  (60)  can  now  be  solved.  The  initial  conditions  used 
were  identical  to  those  used  in  solving  the  quantum  kine- 
tic model.  To  obtain  a solution  the  constants  and  Cg 
are  first  eliminated  from  Eqs  (58)  and  (60)  using  Eqs  (56) 
and  (57).  Eqs  (58)  and  (60)  must  then  be  solved  numeri- 
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The  two  constants  can 


| 

k 


cally  for  the  values  of  T^  and 
then  be  determined  by  substituting  for  T^  and  T^  in  Eqs 
(56)  and  (57)-  The  cascade  growth  rate,  n^<R*>,  is  calcu- 
lated using 


n1C2 


7.V¥(  .)e-*/T2 


and  compared  with  the  corresponding  quantum  kinetic  rate. 

Comparisons  were  made  under  two  sets  of  conditions. 
First,  Eqs  (59)  and  (60)  were  solved  approximately  by  as- 
suming that  most  of  the  energy  went  into  excitations.  The 
third  term  in  Eq  (58)  and  the  third  and  fifth  terms  in  Eq 
(60)  could  then  be  neglected.  Also,  <«> was  assumed  to  be 
(3/2)T^.  The  equations  were  then  resolved  without  any  ap- 
proximations (Appendix  G).  Results  of  the  comparisons 
will  be  discussed  in  Chapter  IV. 


IV.  Analysis  and  Discussion 


The  results  of  each  of  the  comparisons  will  be  pre- 
sented and  discussed  in  this  chapter.  All  comparisons 
were  made  during  the  early  portion  of  the  electron  cas- 
cade, shortly  after  a steady  state  had  been  established. 

The  first  portion  of  the  chapter  discusses  the  results  of 
the  temperature  model  comparisons.  The  remainer  of  the 
chapter  will  deal  with  the  comparisons  involving  the  im- 
proved diffusion  model  and  the  two  temperature  model. 

The  Temperature  Model 

The  comparisons  between  the  temperature  model  and  the 

quantum  kinetic  model  were  made  at  several  values  of  inci- 

7 ■i?  2 

dent  laser  flux  ranging  from  5 x lO'-IO  W/cm  . The  re- 
sults of  these  comparisons  are  presented  in  figs.  2 - 4- 
and  table  I.  Fig.  2 is  a plot  of  the  cascade  growth  rate 
as  a function  of  flux  for  both  the  quantum  kinetic  and  tem- 
perature models.  Figs.  3 and  4-  are  plots  of  the  electron 

8 

distribution  functions  predicted  by  both  models  at  5 x 10 
W/cm2  and  10^W/cm2  respectively.  Table  I summarizes  the 
apportionment  of  energy  among  the  various  loss  mechanisms 
for  several  values  of  incident  flux. 

A comparison  of  the  growth  rates  predicted  by  each  of 
the  models  can  be  made  from  fig.  2.  At  higher  fluxes, 
those  above  5 x lO^W/cm2,  the  rates  predicted  by  the  tem- 
perature model  agree  rather  well  with  the  quantum  kinetic 
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Elec 


Laser  Flux  (W/cm2) 

Pig.  2 Electron  Growth  Rate 


36 


Q p 

rates.  Below  5 x 107W/cm  the  predictions  of  the  two 
models  begin  to  diverge,  disagreeing  by  more  than  four 

Q O 

orders  of  magnitude  at  5 x 107W/cm  . The  rates  predicted 
by  the  temperature  model  are  much  too  large.  This  sug- 
gests that  at  lower  fluxes  a maxwellian  distribution  func- 
tion permits  too  many  electrons  of  higher  energy  to  exist. 
This  can  be  confirmed  using  fig.  3 which  plots  the  two  dis- 

Q p 

tribution  functions  for  an  incident  flux  of  5 x 10  W/cm  . 

At  energies  above  14-  eV,  the  ionization  energy  of  the 
ground  state,  the  quantum  kinetic  distribution  function 
predicts  an  electron  density  significantly  smaller  than 
that  of  the  temperature  model.  Fewer  electrons  are  there- 
fore available  for  ionizations  and  a lower  growth  rate  re- 
sults. In  the  region  between  10.5  eV,  the  ground  state 
excitation  threshold,  and  14-  eV  the  situation  is  reversed. 
The  electron  density  of  the  quantum  kinetic  model  is  now 
greater  than  that  of  the  temperature  model.  Thus,  not 
only  does  the  temperature  model  predict  an  excessively 
high  growth  rate,  but  it  also  permits  fewer  excitations 
than  the  quantum  kinetic  model.  Finally,  below  approxi- 

op 

mately  1.5  x 10  W/cm  the  quantum  kinetic  growth  rate  be- 
comes small  enough  that  many  of  the  neglected  processes 
become  important.  If  these  are  taken  into  account,  the 
quantum  kinetic  growth  rate  would  be  reduced  and  break- 
down might  not  occur.  However,  the  temperature  model  rate 
would  remain  unaffected  and  erroneously  predict  breakdown 
at  fluxes  well  below  the  breakdown  threshold. 
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Electron  Energy  (eV) 

Pig.  3 Electron  Distribution  at  5x10®W/cm2 


Focusing  attention  on  the  high  fluxes,  it  can  be  seen 
that  the  quantum  kinetic  curve  of  fig.  2 begins  to  diverge 
from  the  temperature  model  curve  above  10  W/cm  . Closer 
study  revealed  that  this  divergence  was  artificially  cre- 
ated by  the  conditions  under  which  the  quantum  kinetic  model 
was  solved.  At  sufficiently  high  values  of  incident  flux 
a sizable  number  of  electron  will  possess  energies  in  ex- 
cess of  25  eV;  however,  the  procedure  used  to  obtain  a 
numerical  solution  limits  electrons  to  a maximum  energy  of 
25  eV.  Solution  of  the  quantum  kinetic  model  using  50,  40, 
and  50  1-eV  bins  produced  improved  rates.  As  the  number 
of  bins  was  increased,  the  quantum  kinetic  rates  approach- 
ed those  of  the  temperature  model.  The  effects  of  bin 
size  on  the  quantum  kinetic  rates  was  also  investigated. 

An  increase  to  2 eV  and  3 eV  bins  resulted  in  a signifi- 
cant increase  in  the  growth  rates  while  reducing  the  bin 
size  to  .75  eV  and  .5  eV  had  little  effect  on  the  rates. 
Therefore,  the  choice  of  1 eV  bin  size  appears  to  be 
Justified. 

An  analysis  of  the  energy  gain  anu  loss  mechanisms 
can  be  made  using  table  I.  The  rate  (eV/sec)  at  which  an 
electron  absorbs  energy  from  the  radiation  field  is  given 
by  ft . The  losses  due  to  the  different  mechanisms  are 
given  as  a percentage  of  the  energy  gained.  Losses  which 
result  from  ionization  out  of  the  second  and  third  levels, 
and  excitations  between  these  levels  are  not  included 


since  these  losses  are  10  - 15  orders  of  magnitude  smaller 
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than  the  others  and  need  not  be  considered. 

The  energy  gain  rates  predicted  by  the  two  models 
agree  reasonably  well  over  the  entire  range  of  incident 
fluxes.  This  is  to  be  expected  since  the  classical  micro- 


wave heating  rate  used  in  the  temperature  model  was  de- 
rived directly  from  the  quantum  kinetic  rate  for  net  brem- 
sstrahlung  absorption( Appendix  C).  Although  the  two  rates 
are  independent  of  their  respective  distribution  functions, 
the  small  difference  in  the  rates  can  be  attributed  to  the 


different  average  energies  predicted  by  the  two  models (see 
table  II). 


Table  II 

Average  Electron  Energy  (eV) 


Flux  (W/cm2)  <€> 


TM 

QK 

5 x 107 

1.65 

2.34 

108 

1.80 

2.93 

5 x 108 

2.25 

3.29 

109 

2.45 

3.53 

5 x 109 

3.26 

4.34 

o 

*■ 

o 

(i»l  HW1KT  m 

3.81 

4.86 

5 x 1010 

6.30 

6.90 

1011 

8.84 

9.41 
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At  fluxes  above  5 x 10  W/cm  the  two  models  agree 

rather  well  in  their  predictions  of  energy  lost  to  ground 

I 9 

state  ionization,  R^.  However,  at  fluxes  below  5 x 10' 

p 

W/cm  the  predictions  of  the  two  models  begin  to  diverge, 
with  the  temperature  model  predicting  a significantly  larg- 
er loss  to  ionization  than  the  quantum  kinetic  model.  This 
supports  an  earlier  conclusion,  made  on  the  basis  of  fig. 

3,  that  at  lower  fluxes  the  temperature  model  allows  too 
many  electrons  of  higher  energy  to  exist.  Additional  sup- 
port can  be  found  by  analyzing  eT,  the  energy  required  to 
bring  newly  born  electrons  into  the  distribution  at  the 
average  energy.  At  lower  fluxes  the  larger  temperature 
model  growth  rate  requires  a much  larger  than  does  the 
quantum  kinetic  rate.  This  in  spite  of  the  fact  that  the 
average  quantum  kinetic  energy  is  larger  than  the  average 
temperature  model  energy (see  table  II). 

It  can  be  seen  from  table  I that  relatively  little 
energy  is  lost  to  excitations  between  the  first  and  third 
atomic  levels,  R*^.  This  suggests  that  the  third  atomic 
level  might  be  eliminated  with  negligible  effect  on  the 
growth  rate.  Table  III  summarizes  the  cascade  growth  rates 
predicted  for  both  a two  level  and  three  level  atom.  The 
growth  rates  for  the  two  level  atom  were  obtained  under 
the  same  conditions  as  those  for  the  three  level  atom. 

The  re suit b indicate  an  increase  in  the  growth  rates  as 
might  be  expected.  However,  the  increases  are  small  enough 
that  the  third  level  can  be  neglected  in  most  calculations. 
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Table  III 


1 


Cascade  Growth  Rates  for 
Two  and  Three  Level  Atoms 


Flux 

TM 

QK 

2 

3 

2 

3 

108 

4. 65x1 06 

4-.  4-6x1 06 

2. 94x1 O5 

2.51x10 

5x1 08 

4. 37x1 07 

4.14x1 07 

1 . 66x1 06 

1.44x1 0( 

5x1 09 

l.llxlO9 

1 . 02x1 O9 

7.76x10® 

6 . 86x1  O' 

5x1 010 

2. 98x1 010 

2. 68x1 O10 

2. 94x1 010 

2.70x10' 

1011 

8. 32x1 010 

7. 40x1 O10 

5.71x1010 

5.36x10 

Referring  again  to  table  I,  the  energy  loss  to  exci- 
tations between  the  ground  state  and  second  atomic  level, 
R^2»  can  be  analyzed.  At  the  higher  fluxes  the  two  models 
once  again  predict  similar  losses.  However,  at  lower 
fluxes  the  temperature  model  predicts  a smaller  excitation 
loss  than  the  quantum  kinetic  model.  In  both  models,  ex- 
citations are  the  principal  loss  mechanism  and  account  for 
up  to  90  - 100#  of  the  energy  loss  at  lower  fluxes.  Thus, 
the  lower  excitation  losses  predicted  by  the  temperature 
model  allow  much  more  energy  to  be  available  for  ioniza- 
tions and  increased  growth  rates  result. 

Finally,  the  energy  balance  equation  can  be  rewritten 
as 

R " n^q  “ “T“<R12>n1  (62) 


*3 


where  a two  level  atom  was  assumed  and  R is  the  cascade 
growth  rate.  From  Eq  (62)  it  can  be  seen  that  as  long  as 
the  energy  lost  to  excitations  is  significant  the  growth 
rate  will  depend  strongly  on  the  distribution  function 
through  <R^2>»  Therefore  if  the  growth  rates  predicted 
by  the  two  models  are  to  agree,  their  distribution  func- 
tions  must  be  similar.  At  10  W/cm  fig.  4 shows  the  dis- 
tribution functions  of  the  two  models  to  be  quite  similar 

and  fig.  2 indicates  that  the  predicted  growth  rates  are 

8 2 

almost  identical.  Conversely,  at  5 x 10  W/cm  fig.  3 re- 
veals two  significantly  different  distribution  functions 
and  fig.  2 confirms  an  order  of  magnitude  difference  in 
the  growth  rates. 

The  results  of  the  comparisons  obtained  through  the 
manipulation  of  the  rates  and  energy  levels  will  now  be 
discussed.  Eq  (19) , the  ground  state  ionization  rate  pre- 
dicted by  the  temperature  model,  can  be  rewritten  as 


where  R is  the  cascade  growth  rate  and  a two  level  atom  is 
assumed.  Using  Eq  (63),  the  effect  of  the  manipulations 
on  the  growth  rate  can  be  predicted  in  most  cases.  How- 
ever, since  Eq  (63)  was  obtained  using  the  temperature 
model,  it  cannot  be  expected  to  accurately  predict  quantum 
kinetic  model  results. 

Table  IV  presents  the  results  of  comparisons  made  by 
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manipulating  the  excitation  rate,  R^-,*  ^ (63)  predicts 

that  the  growth  rate  should  increase  as  decreases  and 
decrease  as  R^  increases.  Table  IV  confirms  this  at  all 
values  of  flux  for  both  the  temperature  and  quantum  kinetic 
models.  Eq  (63)  also  accurately  predicts  the  new  values  of 
the  temperature  model  rates.  A closer  comparison  of  the 
rates  in  table  IV  reveals  that  for  .5  Predictions 

O p 

of  the  two  models  agree  at  fluxes  as  low  as  10  W/cm  , while 
for  10  R^  the  rates  of  the  two  models  show  increased  dis- 
agreement. This  is  not  unexpected  and  can  be  explained 
using  energy  considerations.  Since  at  lower  fluxes  almost 
all  the  energy  is  lost  to  excitations,  a decrease  of  50# 
in  the  excitation  cross  section  will  shift  considerable 
energy  into  ionization.  Thus,  the  energy  balance  will  be 
similar  to  that  which  exists  at  higher  fluxes  and  the 
growth  rates,  as  predicted  by  Eq  (62),  will  no  longer  be 
as  strongly  affected  by  small  differences  between  the  two 
distribution  functions.  At  10  R^  the  reverse  will  be 
true.  More  energy  will  be  forced  into  excitations  and  the 
growth  rates  of  the  two  models  will  now  disagree  at  fluxes 
up  to  5 x 10^W/cm^. 

The  results  of  the  comparisons  involving  the  ioniza- 
tion rate,  R*,  are  found  in  table  V.  For  .5  R*  the  growth 
rates  decrease  as  expected  and  the  temperature  model  rates 
are  close  to  the  values  predicted  by  Eq  (63).  For  10  R* 
the  growth  rates  increase,  but  the  temperature  model  rates 
do  not  increase  to  the  extent  predicted  by  Eq  (63).  This 
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happens  because  the  increased  ionization  cross  section 
permits  more  ionizations  at  lower  energies  and  causes  a 

corresponding  decrease  in  the  electron  temperature.  Since 

I I — I/T 

<R  > is  given  approximately  by  RQe  ' , the  decrease  in  T 

will  offset  part  of  the  increase  in  rJ.  A comparison  of 
the  growth  rates  predicted  by  the  two  models  reveals  lit- 
tle change.  At  low  fluxes  the  temperature  model  still 
predicts  a considerably  higher  growth  rate.  This  is  to  be 
expected  since  the  change  in  the  ionization  cross  section 
has  little  effect  on  the  amount  of  energy  going  into  exci- 
tations. At  high  fluxes,  the  agreement  between  the  rates 

predicted  by  the  two  models  is  maintained  with  one  excep- 

I 'll 

tion:  the  quantum  kinetic  growth  rate  for  .5  R/j  at  10 

p 

W/cm  is  too  low.  Reducing  the  ionization  cross  section 
by  50#  increases  the  number  of  electrons  which  possess 
energy  above  14  eV.  To  properly  account  for  this  larger 
number  of  "high  energy  electrons"  the  quantum  kinetic 
equation  must  be  solved  using  an  increased  number  of  bins. 
If  this  is  done  the  quantum  kinetic  growth  rate  for  .5  R^ 
at  10  W/cm  will  approach  the  higher  temperature  model 
rate. 

A change  in  the  ground  state  energy  of  the  atom  will 
affect  the  cascade  growth  rate  through  changes  in  both  the 
ionization  energy  of  the  ground  state  and  the  excitation 
energy  between  the  ground  state  and  the  first  excited 
state.  Table  VI  lists  the  growth  rates  which  result  when 
the  ground  state  energy  is  changed.  Eq  (63)  correctly 
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predicts  the  new  values  of  the  temperature  model  rates. 

An  increase  in  the  ground  state  energy^  = 13)  produces 
a smaller  growth  rate  while  a decrease  in  the  ground  state 
energy (1^  =15)  yields  a larger  rate.  However,  Eq  (63) 
predicts  the  wrong  trend  for  the  quantum  kinetic  rates  at 
5 x 10^  and  lO^W/cm^.  At  these  fluxes  table  VI  shows 
that  an  increase  in  the  ground  state  energy  raises  the 
growth  rate  rather  than  lowering  it  as  predicted  by  Eq 
(63).  Similarly,  a decrease  in  the  ground  state  energy 
produces  a lower  growth  rate.  At  these  fluxes,  a signi- 
ficant amount  of  energy  is  lost  to  ionization  and  the  as- 
sumptions which  were  made  in  deriving  Eq  (63)  lose  their 
validity.  The  change  in  the  growth  rate  becomes  a strong- 
er function  of  1^  and  can  be  better  predicted  from  Eq 
(62).  At  the  lower  fluxes  Eq  (63)  is  valid  and  predicts 
the  correct  trend  of  the  quantum  kinetic  rates.  Comparing 
the  quantum  kinetic  and  temperature  model  growth  rates,  a 
slight  improvement  in  the  agreement  of  the  two  models  is 
noted  for  1^=15.  Conversely,  setting  1^=13  produces  a 
slightly  greater  disagreement  between  the  rates  of  the  two 
models.  This  can  be  explained  by  assuming  that  for  a con- 
stant £ the  electrons  diffuse  along  the  energy  axis  at  a 
constant  velocity.  Let  tj  be  the  time  required  for  an 
electron  to  attain  sufficient  energy(I^)  to  cause  an  ion- 
ization; then  (1-x/j/I^)tI  is  the  fraction  of  tj  during 
which  excitations  can  occur  without  competition  from 
ionization.  Thus,  a smaller  value  of  (l-X/j/I^tj  will 
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allow  more  energy  to  go  into  ionizations.  This  in  turn 
will  yield  better  agreement  between  the  rates  predicted 
by  the  two  models. 

Last  to  be  discussed  are  the  comparisons  obtained  by 
changing  the  first  excited  state  energy.  This  change  will 
affect  only  the  excitation  energy  between  the  ground  state 
and  the  first  excited  state.  The  results  of  the  compari- 
sons are  presented  in  table  VII.  Eq  (63)  correctly  pre- 
dicts both  the  temperature  model  rates  and  the  trend  of 
the  change  in  the  quantum  kinetic  rates.  An  increase  in 
the  excited  state  energy (12=2. 5)  causes  the  growth  rate  to 
increase  while  a decrease  in  the  excited  state  energy 
(l2=4»5)  yields  a smaller  growth  rate.  A closer  compari- 
son of  the  growth  rates  in  table  VII  reveals  that  for 
12*2.5  the  range  of  fluxes  over  which  the  two  models  pre- 
diet  similar  growth  rates  is  extended  to  5 x 10  W/cm  . 

For  I2=4.5  the  disagreement  between  the  two  models  is  in- 
creased at  lower  fluxes,  but  the  agreement  at  higher 
fluxes  is  not  strongly  affected.  This  behavior  can  be 
explained  using  the  same  reasoning  that  applied  to  changes 
in  the  ground  state  energy.  As  the  excitation  energy  in- 
creases, the  time  during  which  excitations  can  occur  with- 
out competition  from  ionization  decreases.  Therefore,  more 
energy  goes  into  ionizations  leading  to  better  agreement 
between  the  growth  rates  predicted  by  the  two  models. 
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The  results  of  the  comparisons  between  the  quantum 
kinetic  and  temperature  models  clearly  pointed  out  the 
need  for  a simple  model  that  could  predict  cascade  growth 
rates  with  reasonable  accuracy  at  lower  fluxes.  Since  Ref 
9 stated  that  the  results  of  the  improved  diffusion  model 
were  valid  over  a wide  range  of  frequencies  and  conditions, 
it  was  hoped  that  this  model  would  be  able  to  predict  the 
electron  growth  rate  over  a wider  range  of  fluxes  than  the 
temperature  model.  Care  was  taken  to  insure  that  the  dif- 
fusion model  was  solved  under  the  same  conditions  as  the 
quantum  kinetic  model.  The  quantity  v*  was  found  to  be 
independent  of  the  exact  form  of  the  distribution  func- 
tion used  over  the  interval  from  xD  to  I . The  growth 
rates  predicted  by  the  diffusion  model  are  given  along 
with  the  corresponding  quantum  kinetic  rates  in  table  VIII. 

Table  VIII 

* 

Diffusion  Model  Growth  Rates 

Flux  QK  * DM 

5 x 108  1.44  x 106  3.90  x 104 

5 x 109  6.86  x 108  2.17  x 10® 

5 x 10™  2.70  x 1010  1.28  x 10™ 


At  5 x 109  and  5 x 10™W/cm^  the  predicted  diffusion  rates 
come  close  to  the  quantum  kinetic  rates;  however,  they  do 


not  appear  to  be  quite  as  accurate  as  the  temperature  model 

Q p 

growth  rates.  At  5 x 10  W/cm  the  rate  predicted  by  the 
diffusion  model  is  two  orders  of  magnitude  less  than  the 
corresponding  quantum  kinetic  rate.  The  predicted  growth 
rate  showed  no  improvement  over  the  temperature  model  rate 
at  the  same  flux.  Finally,  the  failure  of  the  improved 
diffusion  model  to  predict  the  cascade  growth  rate  at  low 
fluxes  appeared  to  be  the  result  of  its  inability  to  accu- 
rately calculate  the  excitation  rate.  Since  this  was  the 
same  problem  experienced  with  the  temperature  model,  con- 
tinued study  of  the  improved  diffusion  model  did  not  ap- 
pear justified. 


The  Two-Temperature  Model 

The  results  of  the  comprisons  between  the  quantum 
kinetic  and  temperature  models  showed  that  at  low  fluxes 
the  actual  electron  distribution  function  could  not  be 


adequately  approximated  by  a maxwellian.  Referring  back 
to  fig.  3,  it  is  seen  that  a maxwellian  allows  too  many 
electrons  to  exist  above  the  ionization  energy  and  too  few 
between  the  excitation  and  ionization  energies.  A close 
study  of  fig.  3 and  similar  plots  at  other  low  values  of 
incident  flux  indicated  that  the  actual  distribution  func- 
tion could  possibly  be  represented  by  two  straight  line 
segments.  These  line  segments,  each  having  a different 
slope,  would  intersect  at  the  excitation  energy  (10.5  eV) 
and  form  a two-temperature  distribution  function.  This 
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forms  the  basis  for  the  two-temperature  model  described 
in  chapter  III. 


i 

l 

I 


jd. 


Since  results  obtained  earlier  in  this  chapter  indi- 
cate that  at  low  fluxes  most  of  the  energy  is  lost  to  ex- 
citations it  was  expected  that  an  approximate  solution  to 
the  two-temperature  model  would  give  results  accurate 
enough  to  yield  a valid  comparison.  The  growth  rates  re- 
sulting from  the  approximate  solution  of  the  two-tempera- 
ture model  are  found  in  table  IX.  A comparison  of  the 

Table  IX 

Two-Temperature  Model  Growth  Rates 


Flux 

QK 

2T 

5 x 107 

1.08  x 

102 

1.70  X 

10“ 

5 x 108 

1.4-4-  x 

106 

4.57  X 

106 

5 x 109 

6.86  x 

108 

7.19  x 

109 

5 x 1010 

2.70  x 

1010 

2.23  x 

101 

quantum  kinetic  and  two-temperature  rates  reveals  that 
only  at  one  flux,  5 x 108W/cm2,  does  the  predicted  two- 
temperature  rate  approach  the  quantum  kinetic  rate.  At 
all  other  fluxes,  including  high  fluxes,  the  predicted 
two-temperature  rates  are  not  even  within  an  order  of 
magnitude  of  the  quantum  kinetic  rates.  The  assumptions 
associated  with  the  solution  were  reexamined  and  the  as- 
sumption that  <«>  * (3/2)T^  was  found  to  be  in  error. 
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Since  T^>T2»  this  assumption  resulted  in  an  overestimate 
of  the  average  energy. 

The  two-temperature  model  was  then  solved  exactly  to 
reduce  the  possibility  of  further  error.  The  resulting 
growth  rates  were,  at  every  value  of  incident  flux,  much 
worse  than  the  rates  of  table  IX.  Agreement  with  the  quan- 
tum kinetic  model  growth  rates  did  not  exist  nor  did  the 
rates  approach  those  of  the  temperature  model  at  high 

fluxes.  Thus,  it  appears  that  the  two-temperature  model, 

< 

which  was  based  on  approximating  the  shape  of  the  distri- 
bution function,  is  not  valid  and  the  possibility  that  the 
electron  distribution  is  a function  of  an  additional  vari- 
able, such  as  incident  laser  flux,  should  be  considered. 
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V.  Conclusion 


Prom  the  results  presented  in  Chapter  IV  several  con- 
clusions can  be  drawn  concerning  the  validity  of  the  vari- 
ous models  assumed. 

Under  the  conditions  assumed  in  Chapter  III  the  tem- 
perature model  adequately  predicts  breakdown  phenomena  for 

q p 

incident  laser  fluxes  above  5 x 1CrW/cm  . At  these  fluxes, 
the  distribution  function  can  be  assumed  to  be  maxwellian. 
In  addition,  at  higher  fluxes,  the  proportion  of  energy 
lost  to  ionizations  is  greater  and  the  growth  rate  is  not 
as  strongly  dependent  on  the  shape  of  the  distribution 

Q p 

function.  Below  5 x 1CrW/cm  the  temperature  model  is  not 
reliable.  The  actual  electron  distribution  is  no  longer 
adequately  approximated  by  a maxwellian  and  the  tempera- 
ture model  growth  rates  are  too  large.  Finally,  a two 
level  atom  can  be  assumed  without  significantly  affecting 
the  accuracy  of  the  temperature  model. 

The  results  of  Chapter  IV  indicate  that  the  range 
over  which  the  quantum  kinetic  and  temperature  model  growth 
rates  agree  can  be  extended.  A change  in  the  excitation 
rate,  ground  state  energy,  and  first  excited  state  energy 
has  a significant  effect  on  the  reliability  of  the  temper- 
ature model.  The  low  flux  accuracy  of  the  temperature 
model  growth  rates  will  improve  as  the  excitation  cross 
section  decreases,  the  ground  state  energy  decreases,  and 
the  first  excited  state  energy  increases.  Therefore,  the 
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temperature  model  will  be  most  reliable  when  applied  under 
the  conditions  of  a low  excitation  cross  section  and  a 
ratio  of  excitation  energy  to  ionization  energy  near  one. 

Attempts  to  develop  a simple  model  which  would  be 
valid  at  low  fluxes  proved  unsucessful.  The  improved  dif- 
fusion model  possessed  sufficient  accuracy  at  high  fluxes, 
but  in  the  limit  of  low  fluxes  it  was  no  better  than  the 
temperature  model.  The  two-temperature  model  could  not 
predict  reasonable  growth  rates  at  any  flux;  therefore,  it 
must  be  concluded  that  the  assumption  of  a two-temperature 
distribution  function  is  not  valid. 

Finally,  it  appears  that  the  key  to  developing  a sim- 
ple breakdown  model  valid  at  low  fluxes  lies  in  a success- 
ful approximation  of  the  electron  distribution  function. 
Therefore,  any  further  investigation  of  breakdown  at  low 
fluxes  should  be  directed  toward  the  development  of  a 
distribution  function  which  will  adequately  approximate 
the  actual  electron  distribution  function.  One  possible 
approach  might  be  a modification  of  the  temperature  model 
distribution  function  to  include  an  additional  term  which 
would  depend  on  both  the  electron  heating  rate  and  the 
electron  energy.  At  high  fluxes  this  term  would  have  lit- 
tle effect  on  the  distribution  function;  however,  at  lower 
fluxes  it  would  modify  the  distribution  function  suffi- 
ciently to  permit  an  accurate  prediction  of  the  cascade 
growth  rate. 
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Appendix  A 


Cross  Sections 

A classical  approach  can  be  used  to  determine  an  ana- 
lytic formula  for  the  ionization  cross  section(Ref  5:392). 
If  an  electron  of  energy  c passes  near  suaother  electron, 
the  differential  cross  section  for  energy  transfer  to  the 
target  electron  in  a range  A«  to  A«  + d(A« ) is  given  by 


d»  - (64) 

(4W«0)2€  (A«)2 

Now,  integrating  this  expression  from  I to  « yields  the 
ionization  cross  section 


*(«) 


ire 


(4"«0)  « 


(65) 


where  I,  the  ionization  energy,  is  the  minimum  energy  that 
can  be  transferred.  Thus,  the  cross  section  for  ioniza- 
tion from  the  jth  level  can  be  written  as 


9 


(4w<)2I, 


(66) 


where  I . represents  the  ionization  energy  of  the  jth  level 

U 

and  ionization  is  assumed  to  occur  each  time  an  energy 
greater  than  1^  is  transferred  to  an  atom. 

A simplified  expression  for  the  excitation  cross 
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a 


section  is  given  by(Ref  10:40) 


Appendix  B 


Rates 


The  rate  for  the  occurrence  of  ionization  or  excita- 
tion will  depend  on  the  cross  section  for  the  interaction 
and  the  relative  velocity  between  the  particles  involved. 
Since  the  velocity  of  the  electrons  will  be  much  greater 
than  that  of  the  atom,  the  velocity  of  the  atom  can  be 
neglected  so  that  the  rate  becomes 

.CO.JE.CO  (33) 

e 


where  < is  the  energy  of  the  electron  and  me  its  mass. 

If  the  cross  section  for  ionization  is  substituted 
into  Eq  (33)  the  rates  of  ionization  used  in  the  quantum 
kinetic  model  are  obtained 


»?<•> 


€ 


we 


e (4tr<o)2lj 


,-3/2 


(«  - I,) 


(34) 


To  obtain  the  ionization  rate,  R*(T),  used  in  the  temper- 
ature model,  Eq  (34)  must  be  averaged  over  a maxwellian 
energy  distribution 


Rj(T)  - | L 7«1/2Rj(«)e"  ^Td«  (70) 

which  upon  integration  becomes 
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Appendix  C 

Derivation  of  the  Temperature  Model 

Using  a method  similar  to  that  of  Ref  11  the  temper- 
ature model  can  be  derived  from  the  quantum  kinetic  equa- 
tion for  the  time  rate  of  change  of  the  electron  distri- 
bution function,  Eq  (23) , where  f(*)  is  assumed  maxwellian. 
For  small  hw/«  , the  bremsstrahlung  absorption  and  emission 
terms  can  be  expanded  in  a Taylor  series  and  in  the  limit 
as  tiuyfc  -*•  0 Eq  (23)  becomes(Ref  2:77^,776) 

■ h , f 

e i f j i 

+ i,f  iniEiJ(,+  Xid)f(*+  ‘ J nJBJCOf(° 

+ rn..R*(«  + Ij)f(«+  Ij) 

♦ 1(0  £ n . /R*(«  )f  (Od«  (71 : 

d d o 3 

I 

where  i is  the  classical  microwave  heating  rate  given  by 

Eq  (9). 

An  equation  for  the  time  rate  of  charge  of  n0  can  be 
Obtained  by  taking  the  zeroth  energy  moment  of  Eq  (71 )• 
Integrating  the  equation  over  all  energy  yields 

“ ne*  nJ<Rd>  W 
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To  obtain  an  equation  for  the  time  rate  of  change  of  <«  > , 
the  average  energy,  the  first  energy  moment  of  Eq  (71)  is 
taken.  Multiplying  the  equation  by  « and  integrating  over 
all  energy  yields 


d<«  > £ ^ t ^ i v 

= - T n .1  .<«  . > 

dt  ne  r o j j 


<«>  dne 


X n.x.  .<R?.> r- XT 

i,J  i 1 ne  dt 


where  the  non-zero  terms  of  the  integration  are 


8f(«)  _ d<«>  . e 

Tt  d<  = ne“3F”  + <€><TF" 


«o  «o 

- /*  iniR?dC* 

+ xid)f(.+  - ne  . £ 


«©  «0 

/«r  n,R*(Of(«)d«  - [*£  n.R^C* 

oj  33  oj  33 


+ Ij)f(«+  Ij)d«  = ne  X 

d 


<R>  being  the  energy  average  of  the  ionization  or  excita- 
tion rate,  and  «>  ■ (3/2)T  for  a maxwellian,  where  T is 
defined  to  be  the  electron  temperature. 

The  equations  for  the  time  rate  of  change  of  the  n^ 
can  be  written  as 


5 


Sr  - * neni  Hi  * dfi<BiJ>nine  + jfi<Ri3>ndne  (75) 

These  three  equations  together  with  Eqs  (4-3)-(44)  form  the 
set  of  equations  which  comprise  the  temperature  model. 
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Appendix  D 


Numerical  Solution  of  the  Quantum  Kinetic  Model 


The  equations  to  he  solved  are  the  quantum  kinetic 
equation  (23)  j and  the  equations  representing  the  changes 
in  the  population  of  the  three  atomic  levels,  Eqs  (39)- 
(41).  The  method  of  solution  employed  is  similar  to  the 
one  used  in  Ref  8. 

Eq  (23) » as  written,  cannot  be  readily  solved,  since 
the  form  of  the  distribution  function  is  unknown.  However, 
if  the  energy  axis  is  divided  into  K discrete  bins  of  width 
fcw  Eq  (23)  is  transformed  into  K equations  of  the  form 

df 

TUT  * naP[am-1fm-1  ' (am+  Vfm  + VlVl] 


+ i,f>iniK^’m+xi/”+xid + 4",',  <76) 


where  fm  represents  the  value  of  the  electron  distribution 
function  at  energy  «m  = (m  - 1/2)hw,  the  midpoint  of  the 
mth  bin.  Rm»  am,  and  bffi  are  the  corresponding  values  for 
the  various  rates  at  energy  «m,  while  I..  and  x^  are  in- 
tegers which  represent  the  number  of  bins  which  correspond 
to  the  various  ionization  and  excitation  energies.  Eq  (76) 
together  with  Eqs  (39)-(41)  can  now  be  solved. 
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The  initial  conditions  are  established  by  specifying 
the  temperature  of  the  gas,  the  photon  flux,  the  number  of 
ground  state  atoms,  the  number  of  electrons  present  as  a 
result  of  non-equilibrium  effects(i.e.  pre-ionization), 
the  number  of  bins,  and  the  bin  width.  The  populations  of 
the  second  and  third  atomic  levels  are  calculated  assuming 
a Boltzmann  distribution 

n.  = n,,  e“x1i/To  (77) 

where  TQ  is  the  gas  temperature  and  the  degeneracy  of  each 
level  is  assumed  to  be  1.  The  total  free  electron  popula- 
tion is  determined  by  adding  the  equilibrium  electron  pop- 
ulation, as  determined  by  the  Saha  equation 


(78) 


to  the  population  resulting  from  non-equilibrium  effects. 
The  electron  population  of  each  bin  is  then  calculated  by 
integrating  a Maxwell-Boltzmann  distribution  function  over 
the  width  of  the  bin 


yW/J  V2 


(79) 


Each  bin  is  assumed  to  have  a minimum  initial  population 
of  one  electron. 

The  rates  for  excitation  and  ionization  are  calculated 
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i 


for  all  bins  which  possess  an  energy,  em,  greater  than  the 
threshold  energy  of  the  corresponding  rate.  They  are  cal- 
culated by  substituting  the  value  of  into  the  appropri- 
ate rates  found  in  table  X.  The  coefficients  of  brems- 
strahlung  absorption  are  calculated  for  each  bin  using 


m 


2e2 

‘2(  cm+  *«)' 

p 

3m  c € w 
e o 

me 

1/2  (cm+W2) 


lu 


m 


(80) 


where  <rm  is  a constant  given  by  Eq  (69).  The  coefficient 
for  bremsstrahlung  emission  is  calculated  for  each  bin  by 
microreversibility 


(81) 


The  emission  and  absorption  rates  are  then  obtained  by 
multiplying  the  emission  and  absorption  coefficients  by 
the  scaled  photon  flux.  The  scaling  of  the  incident  flux 
allows  the  choice  of  a convenient  bin  width,  hw,  which 
does  not  correspond  to  the  actual  photon  energy,  of 

the  incident  radiation.  The  relationship  between  the  ac- 
tual and  scaled  fluxes  is  given  by 

F * — § — 5 h«  (82) 

where  S is  in  W/cm2  and  F is  in  #/(cm2-sec).  The  scaling 
of  Eq  (82)  yields  a correct  net  absorption  rate  for  hw/i 
small.  To  show -this,  the  net  absorption  rate  is  expanded 
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in  a Taylor  series  which  in  the  limit  as  Hw/c-^O  yields 
(Ref  2:776) 


(83) 


where  i is  the  net  bremsstrahlung  absorption  rate  and  S 

p 

is  the  laser  flux  in  W/cm  . 

With  the  initial  conditions  established  Eqs  (39)-(4l) 
and  (76)  are  solved  in  time  using  an  algorithm,  developed 
by  Gear(Ref  12,13) > for  the  solution  of  stiff  systems  of 
first  order  ordinary  differential  equations.  This  scheme 
employs  multistep  predictor-corrector  methods  of  up  to 
sixth  order.  The  error  for  each  step  is  defined  by 


€ = 


(At) 


k+1  (k+1) 
xi 


(k+1  )a*J 


max 


)l/2 


(84) 


where  the  xi's  are  the  solutions  to  the  set  of  equations, 
k is  the  order,  At  is  the  time  step,  x£k+/1^  is  the  (k+1) 
derivative  of  x^,  and  x?^  is  the  maximum  value  of  the  x^ 
encountered.  A step  is  accepted  and  the  time  advanced 
when  the  error  is  less  than  .01  . The  solution  of  the 
equations  is  terminated  once  a specified  value  of  the  free 
electron  density  is  achieved. 
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Appendix  E 


Numerical  Solution  of  the  Temperature  Model 


The  equations  to  be  solved  are  the  equations  for  the 
change  in  population  of  the  three  atomic  levels,  Eqs  (39)- 
(41),  and  the  equations  for  the  zeroth  and  first  energy- 
moments  of  the  Boltzmann  equation,  Eqs  (43)  and  (44).  Eq 
(44)  becomes,  upon  setting  <«>  = (3/2)T, 


(85) 


The  initial  conditions  are  established  in  a manner  quite 
similar  to  that  used  in  the  quantum  kinetic  solution.  The 
gas  temperature,  laser  flux,  density  of  ground  state  atoms, 
and  density  of  non-equilibrium  electrons  are  specified. 

The  populations  of  the  remaining  atomic  levels  are  calcft- 

lated  using  Eq  (77)*  The  total  free  electron  density  is 

- 

determined  by  adding  the  result  of  Eq  (78)  to  the  non- 
equilibrium density.  The  heating  rate  is  calculated  using 
B|  (9).  Finally,  the  rates  as  a function  of  temperature 

ape  calculated  from 


r 
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which  upon  integration  becomes 


RJ(T)  = (C/T5/2)[Te"I/T  + I E±(-I/T)]  (87) 

where  I now  represents  either  the  ionization  or  excitation 
energy,  C is  a constant,  and  E^(-I/T)  is  the  integral  de- 
fined in  Eq  (37)  and  is  evaluated  using  algorithms  found 
in  the  Handbook  of  Mathematical  Functions  (Ref  14:231). 

The  solution  to  the  set  of  equations,  (39)-(4l),  (43) 
and  (83),  is  obtained  by  using  a second  order  predictor- 
corrector  method.  The  equations  are  all  of  the  form 

*i  = ^(Xi>Xj  » * * *Xg)  (88) 

The  various  i^fs  are  first  computed  utilizing  Eqs  (39)-(4l) 
(43)  and  (85),  and  the  initial  conditions  at  tQ.  A pre- 
dicted value  of  x^  at  t = tQ  + At  is  then  computed 

xi  = xi  + *iAt  (89) 

The  derivative  is  computed  and  used  to  obtain  a cor- 
rected value  of  x^  at  t^ 


x?  - (*i  + *?)  (90)- 

The  predicted  and  corrected  values  are  compared  and  if  the 
error  as  defined  by 
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is  less  than  .01,  the  solution  is  accepted,  the  time  ad- 
vanced from  t to  t^ , and  a new  iteration  is  begun  with  At 
doubled.  If  the  error  is  greater  than  .01,  the  solution 
is  rejected,  At  is  reduced  by  a half,  and  the  calculation 
repeated.  This  process  is  then  continued  until  the  solu- 
tion at  a particular  time  is  obtained,  or  a specified 
value  of  any  of  the  x^  is  reached. 


» 
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Appendix  F 


Derivation  of  the  Second  Energy  Moment 
of  the  Boltzmann  Equation 

The  solution  of  the  two-temperature  model  requires 
the  use  of  an  additional  equation.  This  equation  can  be 
obtained  by  taking  the  second  energy  moment  of  the  Boltz- 

O 

mann  equation,  Eq  (11).  Multiplying  Eq  (11)  by  e and 

integrating  over  all  energy  yields  an  equation  for  the 

2 

time  rate  of  change  of  < e > 

= 2~<6>  + x2n1<Rx>  + I2n1<RI> 

■v  T x ^ d.n 

- 2xn>j<«Rx>  - 2In/J<€R1>  - ^_>  ^ (59) 

© 

where  the  non-zero  terms  of  the  integration  are 


? 2 8f  am 

J*  It  ^ = 

d<«2>  ^e 

ne  dt  + <C  > dt 

(92) 

0 

CD 

±-  f.s  ff 

n ^ J 9c  . 

d«  = - 2e<«> 

(93) 

e o 

m 

n^  J«2f(«+  x)RX(c+  x)dc  - n^  J«2f  (« )RX(«  )de  = 

O 0 

x^n^r  Rx>  - 2xn1n0<cRx>  (94) 
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* 

n*  / c2f(c+  l)RI(«  + I)d*  - n1  /«2 f(«)RI(€)  dc  = 

o o 

I2n1ne<RI>  - 2In/,ne<€RI>  (95) 

The  < > describes  the  average  of  the  enclosed  quantity  over 
, the  appropriate  distribution  function. 
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Appendix  G 


Numerical  Solution  of  the 
Two-Temperature  Model 


The  equations  to  be  solved  are  the  equations  for  the 
zeroth,  first,  and  second  energy  moments  of  the  Boltzmann 
equation,  Eqs  (13)-(14)  and  (59),  subject  to  the  condi- 
tions of  Eqs  (55)-(57).  Assuming  a rapid  equilibration  of 

p 

<«>  and  <€  > Eqs  (14)  and  (59)  upon  substitution  for 
dne/dt  become 


xn/,<RX>  - (I  + <e>)n/1<RI>  = 0 


(58) 


2“<«>  + x2n„<Rx>  + (I2  - <e2>)n^<RI> 


-2xn^<«Rx>  - 2In^<«RI>  = 0 


(60) 


Before  Eqs  (58)  and  (60)  can  be  solved,  the  averaged 
quantities  in  each  equation  must  be  evaluated  using  the 
two-temperature  distribution  function,  Eq  (55) • The  aver- 
age value  of  the  energy, <€>  , is  given  by 

C„  / .3/V*T1  d.  ♦ C-  / «5/V */T2  d.  - 


C1T^/2[r(5/2)  - r(5/2,x/T1)]  + C2T|/2r(5/2,x/T2)  (96) 
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where  the  gamma  function  is  defined  by (Ref  14:255) 


r(x)  = 7 tx-1  e_t  dt  (97) 

0 

and  the  incomplete  gamma  function  is  given  by(Ref  14:260) 

r(a,x)  = / ta-1  e_t  dt  (98) 

x 

2 

The  value  of  <«  > is  given  by 

C„  / «5'V*T1  d«  , C,  7.5/2e-  *T2  d«  . 
n 0 x 

C^j 17(7/2)  - r(7/2,x/T1)]  + C2T^/2r(7/2,x/T2)  (99) 

Using  Eqs  (56)  and  (38)  the  average  values  of  the  ioniza- 
tion and  excitation  rate  can  be  written  as 

<RI>  = (2.76x10’7)C2[T2e"I/T2  + I E±(-I/T2)]  (100) 

<RX>  * (6.41x10"7)C2[T2e“x/T2  + x E^-x/T^]  (101) 

where  is  defined  by  Eq  (37).  Finally,  the  averages 
<«R*>  and  <«RX>  upon  integration  become 

<«RJ>  - (2.76x10_7)C2T2  e”I/T2  (102) 

<«RX>  - (6.41x10“7)C2t|  e_x/T2  (103) 
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The  constants  and  C2  must  now  be  expressed  in 
terms  of  T^  and  Solving  Eq  (56)  for  yields 

°1  = C2  e_X/T2  ex/T1  (104) 

Integrating  Eq  (57)  and  substituting  for  using  Eq  (104) 
yields 

C2  = i T1^2  e“x/T2  eX//T1  [r(5/2) 

- r^.x/^)]  + T|/2r(3/2,x/T2))_1  (105) 

Eqs  (58)  and  (60)  can  now  be  solved  for  and  T2 
using  an  iterative  procedure(Ref  15:259).  Eqs  (58)  and 
(60)  are  first  written  as 

P(TvT2)  = ~ - xn1<Rx>  - (I  + <€>)n1<RI>  (106) 

e 

G(T1$T2)  = 2~-<«>  + x2n<1<Rx>  + (I2  - <«2>)n1<RI> 

-2xn1<«Rx>  - 2In1<«RI>  (107) 

The  functions  F(T,|,T2)  and  G(T^,  1'2)  are  then  evaluated  at 

three  points  (T^jTg)^  ^T1,T2^b*  8111(1  ^T1’T2^c  which  are 
estimates  of  the  actual  solution  (T^,T2),  to  Eqs  (58)  and 
(60).  An  improved  estimate,  (T,pT2)d,  is  ot,,tained  using 
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i * 1,2  (108) 


The  improved  estimate  is  used  to  replace  one  of  the  ori- 
ginal estimates  and  the  iteration  is  repeated.  This  pro- 
cedure is  continued  until  successive  estimates  of  (T^,T2) 
agree  to  within  .01  eV.  The  constant  C2  is  then  computed 
using  Eq  (105)  and  the  electron  growth  rate  obtained  using 

R * (2.76x10“7)n1C2[T2e”I/T2  + I Ei(-I/T2)]  (109) 

where  n^  is  the  ground  state  population. 
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